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Abstract 

We review the physics of purely leptonic decays of K^, D^, Df, and pseudoscalar 
mesons. The measured decay rates are related to the product of the relevant weak-interaction- 
based CKM matrix element of the constituent quarks and a strong interaction parameter related 
to the overlap of the quark and antiquark wave-functions in the meson, called the decay constant 
fp. The leptonic decay constants for vr^, K^, D^, Df, and mesons can be obtained with 
controlled theoretical uncertainties and high precision from ab initio lattice-QCD simulations. The 
combination of experimental leptonic decay-rate measurements and theoretical decay-constant cal¬ 
culations enables the determination of several elements of the CKM matrix within the standard 
model. These determinations are competitive with those obtained from semileptonic decays, and 
also complementary because they are sensitive to axial-vector (as opposed to vector) quark flavor¬ 
changing currents. They can also be used to test the unitarity of the first and second rows of the 
CKM matrix. Conversely, taking the CKM elements predicted by unitarity, one can infer “ex¬ 
perimental” values for fp that can be compared with theory. These provide tests of lattice-QCD 
methods, provided new-physics contributions to leptonic decays are negligible at the current level 
of precision. This mini-review was prepared for the Particle Data Group’s 2016 edition, updating 
the versions in Refs. [1-3]. 
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I. INTRODUCTION 


Charged mesons formed from a quark and antiquark can decay to a lepton-neutrino pair 
when these objects annihilate via a virtual W boson. Fig. 1 illustrates this process for the 
purely leptonic decay of a meson. 
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FIG. 1. The annihilation process for pure D~^ leptonic decays in the standard model. 


Similar quark-antiquark annihilations via a virtual to the hnal states occur for 
the 7r~^, D'^, and mesons. (Whenever psuedoscalar-meson charges are specified in 

this article, use of the charge-conjugate particles and corresponding decays are also implied.) 
Let P be any of these pseudoscalar mesons. To lowest order, the decay width is 


r(P ^ £p) 



( 1 ) 


Here Mp is the P mass, is the £ mass, Vq^q^ is the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element between the constituent quarks qiq 2 in P, and Gp is the Fermi coupling 
constant. The decay constant fp is proportional to the matrix element of the axial current 
between the one-P-meson state and the vacuum: 


{0\(lll^ll5q2\P{p)) = iPt^fp , ( 2 ) 

and can be thought of as the “wavefunction overlap” of the quark and antiquark. In this 
article we use the convention in which ^ 130 MeV. 

The decay P^ starts with a spin-0 meson, and ends up with a left-handed neutrino or 
right-handed antineutrino. By angular momentum conservation, the £^ must then also be 
left-handed or right-handed, respectively. In the = 0 limit, the decay is forbidden, and 
can only occur as a result of the finite £ mass. This helicity suppression is the origin of the 
mj dependence of the decay width. Radiative corrections are needed when the final charged 
particle is an electron or muon; for the r they are greatly suppressed due to the large lepton 
mass, and hence negligible. 

Measurements of purely leptonic decay branching fractions and lifetimes allow an ex¬ 
perimental determination of the product \Vq^q^\ fp. If the decay constant fp is known to 
sufficient precision from theory, one can obtain the corresponding CKM element within the 
standard model. If, on the other hand, one takes the value of assuming CKM unitar- 

ity, one can infer an “experimental measurement” of the decay constant that can then be 
compared with theory. 

The importance of measuring F(P —)■ £i>) depends on the particle being considered. 
Leptonic decays of charged pseudoscalar mesons occur at tree level within the standard 
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model. Thus one does not expect large new-physics contributions to measurements of r(P —)■ 
(.u) for the lighter mesons P = and these processes in principle provide clean 

standard-model determinations of Vud and 14s. The situation is different for leptonic decays 
of charm and bottom mesons. The presence of new heavy particles such as charged Higgs 
bosons or leptoquarks could lead to observable effects in r(P —)■ Iv) for P = Pj), [4-8]. 
Thus the determination of |14fe| from P+ —)■ rz/ decay, in particular, should be considered a 
probe of new physics. More generally, the ratio of leptonic decays to tv over ^v hnal states 
probes lepton universality [4, 9]. 

The determinations of CKM elements from leptonic decays of charged pseudoscalar 
mesons provide complementary information to those from other decay processes. The de¬ 
cay P ^ iv proceeds in the standard model via the axial-vector current ?i 7 ^ 75 (? 2 , whereas 
semileptonic pseudoscalar meson decays Pi —)■ P 2 iv proceed via the vector current 
Thus the comparison of determinations of | 14 i<j 2 l from leptonic and semileptonic decays tests 
the V — A structure of the standard-model electroweak charged-current interaction. More 
generally, a small right-handed admixture to the standard-model weak current would lead to 
discrepancies between \Vq^q^\ obtained from leptonic pseudoscalar-meson decays, exclusive 
semileptonic pseudoscalar-meson decays, exclusive semileptonic baryon decays, and inclusive 
semileptonic decays [10, 11]. 

Both measurements of the decay rates r(P —)■ (,v) and theoretical calculations of the 
decay constants fp for P = vr’*', P'"*', from numerical lattice-QCD simulations are now 
quite precise. As a result, the elements of the first row of the CKM matrix |14d| and |14s| can 
be obtained to sub-percent precision from 7r+ —)■ iv and —)■ ^v, where the limiting error 

is from theory. The elements of the second row of the CKM matrix 1 14 ^( 5 )! can be obtained 
from leptonic decays of charged pseudoscalar mesons to few-percent precision, where here 
the limiting error is from experiment. These enable stringent tests of the unitarity of the 
first and second rows of the CKM matrix. 

This review is organized as follows. Because the experimental and theoretical issues asso¬ 
ciated with measurements of pions and kaons, charmed mesons, and bottom mesons differ, 
we discuss each one separately. We begin with the pion and kaon system in Sec. II. First, 
in Sec. IIA we review current measurements of the experimental decay rates. We provide 
tables of branching-ratio measurements and determinations of the product \Vud{s)\fiT+(K+)) 
as well as average values for these quantities including correlations and other effects needed 
to combine results. Then, in Sec. IIB we summarize the status of theoretical calculations of 
the decay constants. We provide tables of recent lattice-QCD results for , fx+, and their 
ratio from simulations including dynamical u, d, s, and (in some cases c) quarks, and present 
averages for each of these quantities including correlations and strong SU(2)-isospin correc¬ 
tions as needed. We note that, for the leptonic decay constants in Secs. IIB, IIIB, and IV B, 
when available we use preliminary averages from the Flavor Lattice Averaging Group [12, 13] 
that update the determinations in Ref. [14] to include results that have appeared since their 
most recent review, which dates from 2013. We next discuss the charmed meson system in 
Sec. Ill, again reviewing current experimental rate measurements in Sec. Ill A and theoret¬ 
ical decay-constant calculations in Sec. IIIB. Last, we discuss the bottom meson system in 
Sec. IV, following the same organization as the two previous sections. 

After having established the status of both experimental measurements and theoretical 
calculations of leptonic charged pseudoscalar-meson decays, we discuss some implications 
for phenomenology in Sec. V. We combine the average i3(P —)■ £v) with the average fp 
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to obtain the relevant CKM elements from leptonic decays, and then compare them with 
determinations from other processes. We also use the CKM elements obtained from leptonic 
decays to test the unitarity of the hrst and second rows of the CKM matrix. Further, as in 
previous reviews, we combine the experimental B{P —>■ f'z/)s with the associated CKM ele¬ 
ments obtained from CKM unitarity to infer “experimental” values for the decay constants; 
the comparison with theory provides a test of lattice and other QCD approaches assuming 
that new-physics contributions to these processes are not signihcant. 


II. PIONS AND KAONS 


A. Experimental rate measurements 


The leading-order expression for r(P —>■ lu) in Eq. (1) is modihed by radiative corrections 
arising from diagrams involving photons, in some cases with additional quark loops. These 
electroweak and “hadronic” contributions can be combined into an overall factor that multi¬ 
plies the rate in the presence of only the strong interaction (T^^^) as follows (cf. Refs. [15, 16], 
and references therein): 


r(p ^ lu) = 


1 + -Cp 

TT 


( 3 ) 


where Cp differs for P = tt,K. The inclusion of these corrections is numerically important 
given the level of precision achieved on the experimental measurements of the vr^ —)■ and 

—)■ decay widths. The explicit expression for the term in brackets above including all 

known electroweak and hadronic contributions is given in Eq. (114) of Ref. [17]. It includes 
the universal short-distance electroweak correction obtained by Sirlin [18], the universal long¬ 
distance correction for a point-like meson from Kinoshita [19] , and corrections that depend 
on the hadronic structure [20]. We evaluate 6p = {a/7i)Cp using the latest experimentally- 
measured meson and lepton masses and coupling constants from the Particle Data group [3] , 
and taking the low-energy constants (LECs) that parameterize the hadronic contributions 
from Refs. [17, 21, 22]. The hnite non-logarithmic parts of the LECs were estimated within 
the large-Ac approximation assuming that contributions from the lowest-lying resonances 
dominate. We therefore conservatively assign a 100% uncertainty to the LECs, which leads 
to a ±0.9 error in We obtain the following correction factors to the individual charged 

pion and kaon decay widths: 


6^ = 0.0176(21) and (5 k = 0.0107(21). 


( 4 ) 


The error on the ratio of kaon-to-pion leptonic decay widths is under better theoretical 
control because the hadronic contributions from low-energy constants estimated within the 
large-Ac framework cancel at lowest order in the chiral expansion. For the ratio, we use the 
correction factor 

Sk/. = -0.0069(17), (5) 

where we take the estimated error due to higher-order corrections in the chiral expansion 
from Ref. [24]. 


^ This uncertainty on is smaller than the error estimated by Marciano and Sirlin in Ref. [23], which 
predates the calculations of the hadronic-structure contributions in Refs. [17, 20-22]. The hadronic LECs 
incorporate the large short-distance electroweak logarithm discussed in Ref. [23] , and their dependence on 
the chiral renormalization scale cancels the scale-dependence induced by chiral loops, thereby removing 
the dominant scale uncertainty of the Marciano-Sirlin analysis [23] . 
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The sum of branching fractions for vr” —)■ and 7r“ —)■ /i“P 7 is 99.98770(4)% [3]. The 

two modes are difficult to separate experimentally, so we use this sum. Together with the 
lifetime 26.033(5) ns [3] this implies r(7r“ —= 3.8408(7) x 10^ s“^. The right-hand 
side of Eq. (1) is modified by the factor 1.0176 ± 0.0021 mentioned above to include photon 
emission and radiative corrections [23, 25]. The decay rate together with the masses from 
the 2014 PDG review [3] gives 

U-\Vud\ = (127.13 ± 0.02 ± 0.13) MeV, (6) 

where the errors are from the experimental rate measurement and the radiative correction 
factor Stt in Eq. (4), respectively. The uncertainty is dominated by that from theoretical 
estimate of the hadronic structure-dependent radiative corrections, which include next-to- 
leading order contributions of Ci(e^p^^) in chiral perturbation theory [17]. 

The data on decays have been updated recently through a global fit to branching ra¬ 
tios and lifetime measurements [26]: B{K~ —)■ p“p[ 7 ]) = 63.58(11)% and tk± = 12.384(15) 
ns. The improvement in the branching ratio is primarily due to a new measurement of 
B{K^ —>■ TT^TT+TT”) from KLOE-2 [27], which is correlated with B{K^ 2 ) through the con¬ 
straint that the sum of individual branching ratios must equal unity. The sum of branching 
fractions for K~ fi~9 and K~ —)■ and the lifetime imply r{K~ —>■ fi~9['y]) = 

5.134(11) X 10^ s“^. Again taking the 2014 PDG masses [3], this decay rate implies 

fK+\Vus\ = (35.09 ±0.04 ±0.04) MeV, (7) 


where the errors are from the experimental rate measurement and the radiative correction 
factor 6k, respectively. 

Short-distance radiative corrections cancel in the ratio of pion-to-kaon decay rates [28]: 

where Sk/tt is given in Eq. (5). The left-hand side of Eq. (8) is 1.3367(28), yielding 

|Ks|/jc- ^ Q 27599 ± 0.00029 ± 0.00024, (9) 

\yud\jw- 

where the first uncertainty is due to the branching fractions and the second is due to 5k/tz- 
Here the estimated error on the hadronic structure-dependent radiative corrections is com¬ 
mensurate with the experimental error. 

In summary, the main experimental results pertaining to charged pion and kaon leptonic 
decays are 


\Vud\fn- = (127.13 ± 0.02 ± 0.13) MeV, (10) 

\Vus\fK+ = (35.09 ± 0.04 ± 0.04) MeV , (11) 

|K.|/ic+ ^ q_27599 ^ 0.00029 ± 0.00024 , (12) 

\yud\jw- 

where the errors are from the experimental uncertainties in the branching fractions and the 
theoretical uncertainties in the radiative correction factors 6p, respectively. 
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TABLE I. Recent lattice-QCD results for f^+, fK+, and their ratio. The upper and lower panels 
show (2 + 1 + l)-flavor and (2 + l)-flavor determinations, respectively. When two errors are shown, 
they are statistical and systematic, respectively. Results for and /k in the isospin-symmetric 
limit mu = m^, are noted with an they are corrected for isospin breaking via Eqs. (13)-(15) 
before computing the averages. Unpublished results noted with a or “!■” are not included in 
the averages. 


Reference 

Nf /.+ (MeV) 

/iy+(MeV) 

Ik+Z fTr+ 

ETM 14 [31]§ 

Fermilab/MILC 14 [32] § 
HPQCD 13 [33]§ 

2+1+1 

2+1+1 

2+1+1 

154.4(1.5)(1.3) 

155.92(13)C^^) 

155.37(20)(28) 

1.184(12)(11) 

1.1956(10)(1?6) 

1.1916(15)(16) 

FLAG 15 average [12, 13]^ 

2+1+1 

155.6(0.4) 

1.193(3) 


RBC/UKQCD 14 [34]*-t 

2+1 

130.19(89) 

155.51(83) 

1.1945(45) 

RBC/UKQCD 12 [35]* 

2+1 

127(3)(3) 

152(3)(2) 

1.199(12)(14) 

Laiho & Van de Water 11 [36]^ 

2+1 

130.53(87)(210) 

156.8(1.0)(1.7) 

1.202(11)(9)(2)(5) 

MILC 10 [37] 

2+1 

129.2(0.4)(1.4) 

156 . 1 ( 4 )(; 6 ) 

1.197(2)(lf) 

BMW 10 [38]* 

2+1 


- 

1.192(7)(6) 

HPQCD/UKQCD 07 [39]* 

2+1 

132(2) 

157(2) 

1.189(2)(7) 


FLAG 15 average [12, 13]^ 

2+1 

130.2(1.4) 

155.9(0.9) 

1.192(5) 

Our average 

Both 

130.2(1.7) 

155.6(0.4) 

1.1928(26) 


§ PDG 2014 value of /,r+ = 130.41(21) MeV used to set absolute lattice scale. 

^ Preliminary numbers shown here may change if further new lattice-QCD calculations are published 
before the deadline for inclusion in the final 2015 FLAG review. 

1 Preprint submitted to Phys. Rev. D. Published RBC/UKQCD 12 results included in Nf = 2-1-1 
average. 

^ Lattice 2011 conference proceedings. 

B. Theoretical decay-constant calcnlations 

Table I presents recent lattice-QCD calcnlations of the charged pion and kaon decay con¬ 
stants and their ratio from simulations with three {Nf = 2 + 1) or four flavors {Nf = 2-I-1-I-1) 
of dynamical quarks. The results have been obtained using several independent sets of gauge- 
held conhgurations, and a variety of lattice fermion actions that are sensitive to different 
systematic uncertainties.^ The lattice-QCD uncertainties on both the individual decay con¬ 
stants and their ratio have now reached sub-percent precision. The SU(3)-breaking ratio 
fK+/fn+ can be obtained with especially small errors because statistical errors associated 
with the Monte Carlo simulations are correlated between the numerator and denominator, 
as are some systematics. The good agreement between these largely independent determi¬ 
nations indicates that the lattice-QCD uncertainties are controlled and that the associated 
error estimates are reliable.^ 

^ See the PDG mini-review on “Lattice Quantum Chromodynamics” [29] for a general review of numer¬ 
ical lattice-QCD simulations. Details on the different methods used in modern lattice-QCD calcula¬ 
tions are provided in Appendix A of the FLAG “Review of lattice results concerning low energy particle 
physics” [14]. 

^ The recent review [30] summarizes the large body of evidence validating the methods employed in modern 
lattice-QCD simulations. 6 












Table I also shows the 2015 preliminary three- and four-flavor averages for the pion and 
kaon decay constants and their ratio from the Flavour Lattice Averaging Group (FLAG) [12, 
13] in the lines labeled “FLAG 15 average.” These preliminary updates of the 2013 FLAG 
averages [14] include only those results from Table I that are published in refereed journals, 
or that are straightforward conference updates of published analyses. In the (2-I-1-I-1)- 
flavor averages, the statistical errors of HPQGD and Fermilab/MILG were conservatively 
treated as 100% correlated because the calculations employ ed some of the same gauge-held 
configurations. The errors have also been increased by the -^/y^/dof to reflect a slight tension 
between the results. There are no four-flavor lattice-QGD results for the pion decay constant 
in Table I because all of the calculations listed use the quantity /^+ to fix the absolute lattice 
scale needed to convert from lattice-spacing units to GeV [31-33]. 

All of the results in Table I were obtained using isospin-symmetric gauge-held configu¬ 
rations, he., the dynamical up and down quarks have the same mass. Most calculations of 
pion and kaon decay constants now include the dominant effect of nondegenerate up- and 
down-quark masses by evaluating the masses of the constituent light (valence) quarks in the 
pion at the physical up- and down-quark masses, respectively, and evaluating the mass of 
the valence light quark in the kaon at the physical m„. Those results obtained with degener¬ 
ate up and down valence quarks are corrected for isospin breaking using chiral perturbation 
theory (yPT) before being averaged. The isospin-breaking corrections at next-to-leading 
order in yPT can be parameterized as [24, 40] 


/tt = fn+ 

fx = fK+ (l — ^SU(2)/2) 
fx _ 1 fx+ 

f-TT \/ 5sU(2) + 1 fTT + 


(13) 

(14) 

(15) 


where the expression for 5su(2) in terms of the quark masses, meson masses, and decay 
constants, is given in Eq. (37) of Ref. [14]. Numerically, values of (5su(2) ~ —0.004 were 
employed by FLAG to obtain the (2-|-l)-flavor averages in Table I, but some direct lattice- 
QGD calculations of (5su(2) give larger values [31, 33, 41] and further studies are needed. 

To obtain the best decay-constant values for comparison with experimental rate mea¬ 
surements and other phenomenological applications, we combine the available (2 -|- 1)- and 
(2 -|- 1 -|- l)-flavor lattice-QGD results, first accounting for the omission of charm sea quarks 
in the three-flavor simulations. The error introduced by omitting charm sea quarks can 
be roughly estimated by expanding the charm-quark determinant in powers of l/rric [42]; 
the resulting leading contribution is of order a* (AQCD/2mc)^ [43]. Taking the MS values 
mdfnc) = 1.275 GeV, Aqcd ~ 340 MeV from FLAG [14], and a{fnc) ~ 0.4, leads to an 
estimate of about 0.7% for the contribution to the decay constants from charm sea quarks. 
The charm sea-quark contribution to ratios of decay constants is expected to be further 
suppressed by the SU(3)-breaking factor (m^ — mdf-^QCD, and hence about 0.2%. 

We can compare these power-counting estimates of charm sea-quark contributions to 
the observed differences between the (2-1-1)- and (2-|-1-1-1)-flavor lattice-QGD averages for 
kaon, Il(s)-meson, and R(s)-decay constants and ratios in Tables I, IV, and VI. Of these, the 
kaon decay constants have been calculated most precisely, and the and three- and four-flavor 
averages for fK+ and fK+ / f-K+ agree within sub-percent errors. Within present uncertainties, 
however, effects of this size in pseudoscalar-meson decay constants cannot be ruled out. 
Therefore, to be conservative, in this review we add in quadrature additional systematic 
errors of 0.7% and 0.2% to all (2-|-l)-flavor decay-constant and decay-constant-ratio averages. 
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respectively, to account for the omission of charm sea quarks. Numerically, this increases 
the errors by at most about 50% for /^+ and less for all other decay constants and ratios, 
indicating that the published (2+l)-flavor lattice-QCD results and uncertainties are reliable. 
Our hnal preferred theoretical values for the charged pion and kaon decay constants are 

Our averages : fT,+ = 130.2(1.7) MeV, fK+ = 155.6(0.4) MeV, = 1.1928(26), (16) 

J-K + 

where f^+ is simply the (2+l)-flavor FLAG average with the error increased by the estimated 
0.7% charm sea-quark contribution. For fK+ and fK+/f-K+i we take a simple weighted 
average of the (2-1-1)- and (2-|-l-|-l)-flavor FLAG values, because they are each obtained 
from a sufficient number of independent calculations that we do not expect there to be 
signihcant correlations. In practice, the addition of the charm sea-quark error has a tiny 
impact on our hnal values in Eq. (16), increasing the uncertainty on /^+ by 0.3 MeV, and 
the central value for fx+ffiT+ by one in the last digit. 


III. CHARMED MESONS 

A. Experimental rate measurements 

Measurements have been made for —)■ p^z/, and —)• Only 

an upper limit has been determined for —)■ Both GLEO-c and BES have made 

measurements of decay using e’''e“ collisions at the ^’(3770) resonant energy where 
D ~pairs are copiously produced. They fully reconstruct one of the H’s, say the D~. 
Gounting the number of these events provides the normalization for the branching fraction 
measurement. They then End a candidate /!■*■, and then form the missing-mass squared, 
= {Eqm — — C^cM — Pfi+) 1 taking into account their knowledge of 

the center-of-mass energy, Eqm, and momentum, pcM, that equals zero in e+e“ collisions. A 
peak at zero inplies the existence of a missing neutrino and hence the /r+z/ decay of the 

D^. GLEO-c does not explicitly identify the muon, so their data consists of a combination 
of and —)■ 7t~^u events. This permits them to do two hts: in one they fit for 

the individual components, and in the other they hx the ratio of r+z///i+z/ events to be that 
given by the standard-model expectation. Thus, the latter measurement should be used for 
standard-model comparisons and the other for new-physics searches. Our average uses the 
hxed ratio value. The measurements are shown in Table 11. 

To extract the value of \Vcd\fD+ we use the well-measured D~^ lifetime of 1.040(7) ps. 
The results include a 1% correction (lowering) of the rate due to the presence of the 
radiative hnal state based on the estimate by Dobrescu and Kronfeld [8]. 

We now discuss the . Measurements of the leptonic decay rate have been made 
by several groups and are listed in Table III [47-53]. We exclude older values obtained 
by normalizing to decay modes that are not well dehned. Many measurements, for 
example, used the 07r+ mode. This decay is a subset of the —)• channel which 

has interferences from other modes populating the K^K~ mass region near the 0, the most 
prominent of which is the /o(980). Thus the extraction of the ehective rate is sensitive 
to the mass resolution of the experiment and the cuts used to define the (j) mass region [54].^ 


^ We have not included the BaBar result for B{D'^ —>■ reported in Ref. [55] because this measurement 

determined the ratio of the leptonic decay rate to the hadronic decay rate r(i2+ —>■ t^v)/T{D'^ —)■ 




To find decays in the signal channels, CLEO, BaBar and Belle rely on fully recon¬ 
structing all the hnal state particles except for neutrinos and using a missing-mass technique 
to infer the existence of the neutrino. CLEO uses —?■ D^D* collisions at 4170 MeV, 
while Babar and Belle use e’''e“ —?• DKmrDl collisions at energies near the T(4S'). CLEO 
does a similar analysis as was done for the above. Babar and Belle do a similar 
calculation by using the reconstructed hadrons, the photon from the D*'^ decay and a de¬ 
tected To get the normalization they do a £t without the p’*' and use the signal 

at the Df mass squared to determine the total Df yield. 

When selecting the —>■ vr+n and —)■ p+P decay modes, CLEO uses both the cal¬ 

culation of the missing-mass and the fact that there should be no extra energy in the 
event beyond that deposited by the measured tagged D~ and the decay products. The 
—)■ mode, however, uses only extra energy. Babar and Belle also use the extra 

energy to discriminate signal from background in their measurements. 

We extract the decay constant times the CKM factor from the measured branching ratios 
using the mass of 1.96830(11) GeV, the mass of 1.77682(16) GeV, and a lifetime of 
0.500(7) ps [3]. GLEO has included the radiative correction of 1% in the rate listed 
in the Table [8] (the rates need not be corrected). Other theoretical calculations show 
that the rate is a factor of 40-100 below the rate for charm [57-66]. As this is a 
small effect we do not attempt to correct the other measurements. The values for fjy+\Vcs\ 
are in good agreement for the two decay modes. Our average value including both the 
and hnal states is 250.9 ± 4.0 MeV. 


B. Theoretical decay-constant calcnlations 

Table IV presents recent theoretical calculations of the charged D~^- and Hs-meson decay 
constants and their ratio. The upper two panels show results from lattice-QGD simulations 
with three {Nf = 2 -|- 1) or four havors {Nf = 2 -|- 1 -|- 1) of dynamical quarks. Although 
there are fewer available results than for the pion and kaon sector, both fij+ and have 
been obtained using multiple sets of gauge-held conhgurations with diherent lattice fermion 
actions, providing independent conhrmation. For comparison, the bottom panel of Table IV 
shows non-lattice determinations from QGD sum rules and the light-front quark model; only 
results which include uncertainty estimates are shown. The lattice and non-lattice results 


TABLE II. Experimental results for B{D^ —)• B{D~^ —)• r+n), and \Vcd\fD+- Numbers for 

\^cd\fD+ have been extracted using updated values for masses (see text). Radiative corrections are 
included. Systematic uncertainties arising from the Z)’*' lifetime and mass are included. For the 
average //"’“n number we use the CLEO-c result for -|- r'’'n. 


Experiment 

Mode 

B \Vcd\fD+ (MeV) 

CLEO-c [44, 45] 


(3.93 ± 0.35 ± 0.09) x 10“^ 47.07 ± 2.10 ± 0.57 

CLEO-c [44, 45] p+n+r+n 

(3.82 ± 0.32 ± 0.09) x 10"^ 46.41 ± 1.94 ± 0.57 

BES [46] 

p+n 

(3.71 ± 0.19 ± 0.06) X 10"^ 45.73 ± 1.17 ± 0.38 

Our average 

Lines 2-|-3 

(3.74 ± 0.17) X 10"^ 45.91 ±1.05 

CLEO-c [47, 48] 

r+n 

< 1.2 X 10-3 
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TABLE III. Experimental results for B{D^ —> and \Vcs\fjj+- Numbers 

for \Vcs\fD+ have been extracted using updated values for masses (see text). The systematic 
uncertainty for correlated error on the Df lifetime is included. The mass uncertainties are also 
common, but negligible. Common systematic errors in each experiment have been taken into 
account in the averages. 


Experiment 

Mode 

R(%) 

\Vcs\fj,+ (MeV) 

CLEO-c [47, 48] p+u 

(0.565 ± 0.045 ± 0.017) 250.8 ± 10.0 ± 4.2 

BaBar^ [53] 


(0.602 ± 0.038 ± 0.034) 

258.9 ± 8.2 ±7.5 

Belle [49] 


(0.531 ±0.028 ±0.020) 

243.1 ± 6.4 ±4.9 

Our average 


(0.556 ±0.024) 

248.8 ± 5.8 


CLEO-c [47, 48] 

r'*“zz (tt+E) 

(6.42 ±0.81 ±0.18) 

270.8 ± 17.1 ±4.2 

CLEO-c [50] 

r"'“zz (p'^u) 

(5.52 ±0.57 ±0.21) 

251.1 ±13.0 ±5.1 

CLEO-c [51, 52] 

T'^i' {e^oT') 

(5.30 ±0.47 ±0.22) 

246.1 ±10.9 ±5.4 

BaBar [53] 

T~^l' {pL^)vi') 

(5.00 ±0.35 ±0.49) 

239.0 ± 8.4 ± 11.9 

Belle [49] 


(6.04 ± 0.43l°:^®) 

262.7 ± 

Belle [49] 

T~^v {e'^vV) 

(5.37 ± 0.33l[J|f) 

247.7 ±7.6^^;^ 

Belle [49] 


(5.86 ± 0.37+°|4) 

258.7 ± 8.21^3^2 


Our average 


(5.56 ±0.22) 

252.1 ±5.2 

Our average 



250.9 ±4.0 


We do not use a previous unpublished BaBar result from a subsample of data that uses a different 
technique for obtaining the branching fraction normalization [56]. 


agree, but the uncertainties on D^^^-meson decay constants from lattice QCD have now 
reached signihcantly greater precision than those from other approaches. 

The lattice-QCD results in Table IV were all obtained using isospin-symmetric gauge-held 
conhgurations. The two calculations by the Fermilab Lattice and MILC Collaborations [32, 
69], however, include the dominant strong isospin-breaking contribution by evaluating the 
mass of the valence light quark in the T)+-meson decay constant at the physical down-quark 
mass. Reference [32] provides a determination of the size of this correction, 

/D+-/D = 0.47(l)Cf)MeV. (17) 

where fo is the value of the D-meson decay constant evaluated at the average up-down 
quark mass. Equation (17) implies that the correction to the SU(3)/-breaking ratio is 

^ =-0.0026, (18) 

JD+ Id 

taking the central values for fjj+ and from the same work. Because the errors on the 
calculations listed in Table IV that neglect isospin breaking are still about 5-8 x larger 
than the sizes of the shifts in Eqs. (17)-(18), we do not correct any results a posteriori for 
this effect in the current review. Nevertheless, we strongly encourage future lattice-QCD 
publications to present results for both the D^- and D°-meson decay constants. Including 
the effect of isospin breaking will be essential once lattice-QCD calculations of fo and fujln 
reach the level of precision in Eqs. (17)-(18). 
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TABLE IV. Recent theoretical determinations of /£>+, fos^ ^ind their ratio. The upper panels 
show results from lattice-QCD simulations with (2 + 1 + 1) and (2 + 1) dynamical quark flavors, 
respectively. Statistical and systematic errors are quoted separately. Lattice-QCD results for Jd 
and fos/fo in the isospin-symmetric limit ruu = Trifi are noted with an The bottom panel 
shows estimates from QCD sum rules (QCD SR) and the light-front quark model (LFQM). These 
are not used to obtain our preferred decay-constant values. 


Reference 

Method 

Nf 

/D+(MeV) 

/D%MeV) 

JDs/ fD+ 

ETM 14 [31]* 

LQCD 

2+1+1 

207.4(3.7)(0.9) 

247.2(3.9)(1.4) 

1.192(19)(11) 

Fermilab/MILC 14 [32] 

LQCD 

2+1+1 

212.6(0.4)(1}:0) 

249.0(0.3)(;ji) 

1.1712(10)(11) 

Average 

LQCD 

2+1+1 

212.2(1.5) 

248.8(1.3) 

1.172(3) 

XQCD 14 [67]* 

LQCD 

2+1 

- 

254(2)(4) 

- 

HPQCD 12 [68]* 

LQCD 

2+1 

208.3(1.0)(3.3) 

- 

1.187(4)(12) 

Fermilab/MILC 11 [69] 

LQCD 

2+1 

218.9(9.2)(6.6) 

260.1(8.9)(6.1) 

1.188(14)(21) 

HPQCD 10 [70]* 

LQCD 

2+1 

- 

248.0(1.4)(2.1) 

- 

Average 

LQCD 

2+1 

209.2(3.3) 

249.8(2.3) 

1.187(12) 

Our average 

LQCD 

Both 

211.9(1.1) 

249.0(1.2) 

1.173(3) 

Wang 15 [71]^ 

Gelhausen 13 [72] 
Narison 12 [73] 

QCD SR 
QCD SR 
QCD SR 


208(10) 

201 ra 

204(6) 

240(10) 

238(1^3) 

246(6) 

1.15(6) 

i-isCT 

1.21(4) 

Lucha 11 [74] 

QCD SR 


206.2(8.9) 

245.3(16.3) 

1.193(26) 

Hwang 09 [75] 

LFQM 


- 

264.5(17.5)’ 

1.29(7) 


§ Obtained using results using are also given in the paper. 

^ Obtained by combining PDG value fo = 205.8(8.9) MeV [76] with fojfo from this work. 


We average the lattice-QCD results in Table IV accounting for possible correlations be¬ 
tween them following the approach established by Laiho et al. [77]. Whenever we have 
reason to believe that a source of uncertainty is correlated between two results, we conserva¬ 
tively take the correlation to be 100% when calculating the average. We then construct the 
correlation matrix for the set of lattice-QCD results using the prescription of Schmelling [78]. 

We hrst separately average the three- and four-flavor results for the charged Dj^-meson 
decay constants and their ratio. There have been no new three-flavor lattice-QCD calcula¬ 
tions of /d+ or /^+//d+ since 2013, so we take the (2+l)-flavor averages from FLAG [14]. 
In this average, the statistical errors were treated as 100% correlated between the results of 
Fermilab/MILC [69] and HPQCD [68] because the calculations employed some of the same 
ensembles of gauge-held conhgurations. For we average the (2+l)-havor results given 
in Table IV, again treating the Fermilab/MILC [69] and HPQCD [70] statistical errors as 
correlated, and taking the xQCD result [67] to be independent. For the (2 + 1 + l)-havor 
D(s)-meson decay constants, we take a simple weighted average of the ETM [31] and Fer¬ 
milab/MILC 14 results [32] in Table IV. We expect them to be independent because the 
calculations use different light-quark and gluon actions and different treatments of the chiral- 
continuum extrapolation. Our separate three- and four-havor averages are listed in the lines 
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labeled “Average” in Table IV, where the errors on the ( 2 +l)-flavor and ( 2 +l+l)-flavor 
fn averages have been rescaled by the factors \/(x^/dof) = 1.1 and \/(x^/dof) = 1.3, 
respectively.^ 

To obtain the single-best values of the H^^-meson decay constants for phenomenology 
applications, we combine the available (2 -|- 1)- and (2 -|- 1 -|- l)-flavor lattice-QCD results, 
which are compatible within the current level of precision. We account for the omission of 
charm sea-quark contributions in the three-flavor calculations by adding to the errors on the 
(2-|-l)-flavor averages in Table IV our power-counting estimates of charm sea-quark errors 
from Sec. IIB. Because the estimated charm sea-quark errors of 0.7% for decay constants 
and 0 . 2 % for decay-constant ratios are less than those on the ( 2 -|-l)-flavor averages, adding 
them in quadrature has a small impact on the total uncertainties. The error increase is at 
most about 25% for and below 10% for both /^+ and fojfD+- Our hnal preferred 
theoretical values for the charged Hj^-meson decay constants are given by the weighted 
average of the entries in the two lines labeled “Average” in Table IV, after including the 
additional charm sea-quark errors in the ( 2 -|-l)-flavor entries: 

Our averages : /b+ = 211.9(1.1) MeV, fn, = 249.0(1.2) MeV , ^ = 1.173(3). (19) 

JD+ 

In practice, the errors on the (2-|-1-|-1)-flavor averages are so much smaller than on the 
(2-1-1) -flavor averages that the combination in Eq. (19) is almost identical to the (2-I-1-I-1)- 
flavor average in Table IV. The most precise result from Fermilab/MILC, in particular, has 
a large weight in the average. 


IV. BOTTOM MESONS 

A. Experimental rate measurements 

The Belle and BaBar collaborations have found evidence for B~ —)■ t~v decay in e’''e“ —)■ 
B~ collisions at the T(4S') energy. The analysis relies on reconstructing a hadronic 
or semi-leptonic B decay tag. Ending a r candidate in the remaining track and photon 
candidates, and examining the extra energy in the event which should be close to zero for a 
real t~ decay to e~ vv or opposite a B~^ tag. While the BaBar results have remained 

unchanged. Belle reanalyzed both samples of their data. The branching fraction using 
hadronic tags changed from 1-79 x lO”"^ [80] to 0.72 +q; 25 ± 0.11 x lO”"^ [81], while 

the corresponding change using semileponic tags was from 3 ® to 1.25±0.28±0.27. 

These changes demonstrate the difficulty of the analysis. The results are listed in Table V. 

There are large backgrounds under the signals in all cases. The systematic errors are 
also quite large. Thus, the significances are not that large. Belle quotes 4.On for their 
combined hadronic and semileptonic tags, while BaBar quotes 3.3(T and 2.3 a, for hadronic 
and semileptonic tags. Greater precision is necessary to determine if any effects beyond the 
Standard Model are present. 

To extract the value of \Vub\fB+ we use the PDG 2014 value of the B~^ lifetime of 
1.638±0.004 ps, and the and B^ masses of 1.77684 and 5.27926 GeV, respectively. 


® After this article was submitted for review, preliminary (2-1-1)- and (2-|-l-|-l)-flavor FLAG averages for 
Id, fos, and fnjfo were presented in Ref. [79] that are identical to our separate averages in Table IV. 
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TABLE V. Experimental results for B{B —)• r n) and \ Vub\fB+- 


Experiment Tag 

B (units of 10 '^) |14b|/s+ (MeV) 

Belle [81] 
Belle [82] 
Belle [82] 

Hadronic 0.72'^q'H ± 0.11 

Semileptonic 1.25 ± 0.28 ± 0.27 
Average 0.91 ± 0.22 

0.72 + 0.09 

BaBar [83] 
BaBar [84] 
BaBar [83] 

Hadronic 

Semileptonic 

Average 

1.830.24 

1.7 ±0.8 ±0.2 

1.79 + 0.48 

1.01 + 0.14 

Our average 


1.06 + 0.20 

0.77 + 0.07 


B. Theoretical decay-constant calcnlations 


Tables VI and VII present theoretical calculations of the and i?s-meson decay 

constants and their ratios. (The decay constants of the neutral B^ and Bg mesons enter the 
rates for the rare leptonic decays B^^g —t ■) The upper two panels show results from 

lattice-QCD simulations with three {Nj = 2 + 1) or four flavors {Nf = 2 + 1 + 1) of dynamical 
quarks. For all decay constants, calculations using different gauge-field configurations, light- 
quark actions, and 6 -quark actions provide independent confirmation. For comparison, the 
bottom panel of Table VI shows non-lattice determinations of the i?(s)-meson decay constants 
which include error estimates. These are consistent with the lattice values, but with much 
larger uncertainties. 


The lattice-QCD results in Tables VI and VII were all obtained using isospin-symmetric 
gauge-field configurations. The most recent calculations of /b+ by the HPQCD, Fermi- 
lab/MILC, and RBC/UKQCD Collaborations [69, 86 , 88 ], however, include the dominant 
effect of nondegenerate up- and down-quark masses by evaluating the decay constant with the 
valence light-quark mass fixed to the physical up-quark mass. HPQCD and RBC/UKQCD 
also calculate fso by fixing the valence light-quark mass equal to the physical down-quark 
mass [ 86 , 88 ]; they find differences between the B~^- and R°-meson decay constants of 
fso ~ fB+ ~ 4 MeV and fBjfB+ ~ Ibs/ fB+ ~ 0.025. Inspection of Tables VI and VII shows 
that these differences are comparable to the error on the HPQCD 12 result for [89], and 
to the errors on the Fermilab/MILC, HPQCD 12, and ETM results for fBs/fs [85, 89], none 
of which account for isospin-breaking. Therefore, to enable comparison with experimental 
measurements, in this review we correct those lattice-QCD results for R-meson decay con¬ 
stants obtained with degenerate up and down valence quarks a posteriori for isospin breaking 
before computing our averages. For the correction factors, we use the differences obtained 
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TABLE VI. Recent theoretical determinations of /b+, Jbs: ^ind their ratio. The upper panels 
show results from lattice-QCD simulations with (2 + 1 + 1) and (2 + 1) dynamical quark flavors, 
respectively. For some of the lattice-QCD results, statistical and systematic errors are quoted 
separately. Lattice-QCD results for fs and fsjfs in the isospin-symmetric limit are 

noted with an they are corrected by the factors in Eqs. (20) and (21), respectively, before 
computing the averages. Preliminary conference results noted with a “i” are not included in the 
averages. The bottom panel shows estimates from QCD sum rules and the light-front quark model, 
which are not used to obtain our preferred decay-constant values. 


Reference 

Method 


/n+(MeV) 

/iJ.(MeV) 

IBs/ fB+ 

ETM 13 [85]*’t 

LQCD 

2 + 1+1 

196(9) 

235(9) 

1.201(25) 

HPQCD 13 [ 86 ] 

LQCD 

2 + 1+1 

184(4) 

224(5) 

1.217(8) 

Average 

LQCD 

2 + 1+1 

184(4) 

224(5) 

1.217(8) 

Aoki 14 [87]*4 

LQCD 

2+1 

218.8(6.5)(30.8) 

263.5(4.8)(36.7) 

1.193(20)(44) 

RBC/UKQCD 14 [ 88 ] 

LQCD 

2+1 

195.6(6.4)(13.3) 

235.4(5.2)(11.1) 

1.223(14)(70) 

HPQCD 12 [89]* 

LQCD 

2+1 

191(1)(8) 

228(3)(10) 

1.188(12)(13) 

HPQCD 12 [89]* 

LQCD 

2+1 

189(3)(3)* 


- 

HPQCD 11 [90] 

LQCD 

2+1 

-■ 

225(3)(3) 

- 

Fermilab/MILC 11 [69] 

LQCD 

2+1 

196.9(5.5)(7.0) 

242.0(5.1)(8.0) 

1.229(13)(23) 

Average 

LQCD 

2+1 

189.9(4.2) 

228.6(3.8) 

1.210(15) 

Our average 

LQCD 

Both 

187.1(4.2) 

227.2(3.4) 

1.215(7) 

Wang 15 [71]^ 

QCD SR 


194(15) 

231(16) 

1.19(10) 

Baker 13 [91] 

QCD SR 


186(14) 

222 ( 12 ) 

1.19(4) 

Lucha 13 [92] 

Gelhausen 13 [72] 
Narison 12 [73] 

QCD SR 
QCD SR 
QCD SR 


192.0(14.6) 

207(^1^) 

206(7) 

228.0(19.8) 

242 (t 17) 
234(5) 

1.184(24) 

1.17(t3) 

1.14(3) 

Hwang 09 [75] 

LFQM 


- 

270.0(42.8)’ 

1.32(8) 


i Lattice 2013 conference proceedings, 
i Obtained with static b quarks (ie. mi, —>■ oo). 

* Obtained by combining from HPQCD 11 with Jbs/ Ib from this work. Approximate statistical 
(systematic) error obtained from quadrature sum of individual statistical (systematic) errors. 

^ Obtained using results using m^°^^ are also given in the paper. 

^ Obtained by combining PDG value Jb = 204(31) MeV [76] with /b,/ fs from this work. 


empirically by HPQCD in Ref. [86]® 


/b+-/ s =-1.9(5) MeV, 

(20) 

hs hs ^ q_qi2 (4) , 

Jb+ Jb 

(21) 

/bo -/s = 1.7(5) MeV, 

(22) 

= -0.011(4). 

/bo Jb 

(23) 


® The correlated uncertainties were provided by HPQCD via private communication. 
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TABLE VII. Recent lattice-QCD determinations of f^o and /bs/Ib^- Results obtained in the 
isospin-symmetric limit = rrid are noted with an while those for the R+-meson are noted 
with an Although the quoted results are identical those in Table VI, they are corrected by 
different factors in Eqs. (20)-(23) before computing the averages. Other labels and descriptions 
are the same as in Table VI. 


Reference 

Method 

Nf 

/so(MeV) 

fBs/fBO 

ETM 13 [85]* *d 

LQCD 

2+1+1 

196(9) 

1.201(25) 

HPQCD 13 [86] 

LQCD 

2+1+1 

188(4) 

1.194(7) 

Average 

LQCD 

2+1+1 

188(4) 

1.194(7) 

Aoki 14 [87]*d 

LQCD 

2+1 

218.8(6.5)(30.8) 

1.193(20)(44) 

RBC/UKQCD 14 [88] 

LQCD 

2+1 

199.5(6.2)(12.6) 

1.197(13)(49) 

HPQCD 12 [89]* 

LQCD 

2+1 

191(1)(8) 

1.188(12)(13) 

HPQCD 12 [89]* 

LQCD 

2+1 

189(3)(3)* 

- 

Fermilab/MILC 11§ [69] 

LQCD 

2+1 

196.9(5.5)(7.0) 

1.229(13)(23) 

Average 

LQCD 

2+1 

193.6(4.2) 

1.187(15) 

Our average 

LQCD 

Both 

190.9(4.1) 

1.192(6) 


^ Lattice 2013 conference proceedings. 

^ Obtained with static b quarks (ie. to;, — >■ oo). 

* Obtained by combining from HPQCD 11 with JeJIb from this work. Approximate statistical 
(systematic) error obtained from quadrature sum of individual statistical (systematic) errors. 


The isospin-breaking correction factors in Eqs. (20)-(23) are well determined because of 
cancellations between correlated errors in the differences. 

We first average the published (2-1-1)-flavor lattice-QCD results for the charged and neu¬ 
tral i?(s)-meson decay constants and their ratios in Tables VI and VII, accounting for possibly 
correlated uncertainties. We treat the statistical errors as correlated between the calcula¬ 
tions of Aoki et al. and RBC/UKQCD because they employ the same gauge-held conhgu- 
rations [87, 88].^ We also treat the statistical errors as correlated between the HPQCD and 
Fermilab/MILC calculations because they analyze an overlapping set of gauge-held conhgu- 
rations [69, 89, 90]. For we include HPQCD’s results from both 2011 [90] and 2012 [89], 
which were obtained using diherent 6-quark actions, but on some of the same gauge-held 
conhgurations. HPQCD 11 and 12 also use the same determination of the absolute lattice 
scale, which is the second-largest source of systematic uncertainty in both calculations. We 
therefore treat the statistical and scale errors as correlated between HPQCD’s (2-|-l)-havor 
Jb^ results. HPQCD also presents two results for in Ref. [89]. The more precise value is 
obtained by combining the ratio fBs/fs from this work with Jb^ from Ref. [90], but an asso¬ 
ciated error budget is not provided. Because this would be needed to estimate correlations 
between the two /s determinations, we include only HPQCD’s more precise (2-|-l)-havor re¬ 
sult for Jb in our average. Our separate three- and four-havor averages for the B~^-, B^-, and 

^ There may be mild correlations between some sub-dominant systematic errors of Aoki et al. and 
RBC/UKQCD, who use the same determinations of the absolute lattice scale and the physical light- 
and strange-quark masses from Ref. [93], and who use the same power-counting estimates for the light- 
quark and gluon discretization errors. The effects of any correlations between these systematics, however, 
would be too small to impact the numerical values of the averages. 
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i?s-nieson decay constants and ratios are listed in the lines labeled “Average” in Tables VI 
and VII, where the error on the (2+l)-flavor average has been rescaled by the factor 
■\/(x^/dof) = 1.2 to account for the tension among results. Our (2+l+l)-flavor “averages” 
are identical to the “HPQCD 13” entries in Tables VI and VII, whcih are the only published 
four-flavor results available. 

To obtain the single-best values of the i?(s)-meson decay constants for phenomenology 
applications, we combine the available (2-1-1)- and (2-|-1-1-1)-flavor lattice-QCD results, which 
are compatible within the current level of precision. Because the four-flavor “average” is 
obtained from only a single result, we do not simply combine the two lines labeled “Average” 
in Tables VI and VII, which would weight the four-flavor result too heavily. Instead, we form 
a single average including the published (2-|-l)-flavor results and the (2-|-l-|-l)-flavor result 
from HPQCD 13. We account for the omission of charm sea-quark contributions in the 
three-flavor calculations by adding to the errors on the (2-|-l)-flavor averages in Tables VI 
and VII our power-counting estimates of charm sea-quark errors from Sec. IIB, taking 
charm sea-quark error to be 100% correlated between the three-flavor results. Because the 
estimated charm sea-quark errors of 0.7% for decay constants and 0.2% for decay-constant 
ratios are much less than those on the (2-1-1)-flavor averages, adding them in quadrature 
has a tiny impact on the total uncertainties. The largest observed change is an 0.3 MeV 
increase on the error from HPQCD 11, and most are negligible. In the combined three- 
and four-flavor average we also consider correlations between the results of HPQCD 12 and 
HPQCD 13 because, although they employ different gauge-held conhgurations, they both 
use NRQCD for the 6-quark action and the bottom-light axial-vector current.® We take 
both the operator-matching and relativistic errors, which are the dominant uncertainties 
in the decay constants, to be correlated between the two calculations. Our hnal preferred 
theoretical values for the charged and neutral -meson decay constants and their ratio 
are 

Our averages : fB+ = 187.1(4.2) MeV, /s, = 227.2(3.4) MeV, ^ = 1.215(7), (24) 

JB+ 

fn 

Jbo = 190.9(4.1) MeV, ^ = 1.192(6). (25) 

JBO 

The errors on f^, /^, and Jb^ after combining the three- and four-havor results are only 
slightly smaller than those of the separate averages due to the correlations assumed. 


V. PHENOMENOLOGICAL IMPLICATIONS 
A. \Vud\i \yus\i and status of first-row unitarity 

Using the average values for / 7 r+|Uud|, /ic+|I4s|, and their ratio from Eqs. (10)-(12) and 
for 7,^+, fK+, and their ratio from Eq. (16), we obtain the following determinations of the 
CKM matrix elements iWdl, |Ws|; and their ratio from leptonic decays within the standard 
model: 

iKrfI = 0.9764(2)(127)(10), |K,| = 0.2255(3)(6)(3), ^ = 0.2314(2)(5)(2), (26) 

I Wdl 

® HPQCD 13 uses a 1-loop radiatively improved 6-quark action, whereas HPQCD 12 uses tree-level action 
coefficients. Both include the same contributions to the currents at one loop, but renormalisation details 
differ. 
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where the errors are from the experimental branching fraction(s), the pseudoscalar decay 
constant (s), and radiative corrections, respectively. These results enable a precise test of the 
unitarity of the first row of the CKM matrix from leptonic decays alone (the contribution 
from |14fe| is negligible). Using the values of \Vud\ and \Vus\ from Eq. (26), we hnd 

iKdP + |K,|2 + |K,|2 - 1 = 0.004(25), (27) 

which is consistent with three-generation unitarity at the sub-percent level. 

The determinations of |14d| and \Vus\ from leptonic decays in Eq. (26) can be com¬ 
pared to those obtained from other processes. The result above for |14d| agrees with 
the determination from superallowed /3-decay, \Vud\ = 0.97417(21) [94], but has an er¬ 
ror more than hfty times larger that is primarily due to the uncertainty in the theoreti¬ 
cal determination of f.,r+- The CKM element |14s| can be determined from semileptonic 
—)■ decay. Here experimental measurements provide a value for the product 

f^^{0)\Vus\, where is the form-factor at zero four-momentum transfer between the 

initial state kaon and the hnal state pion. Taking the most recent experimental determina¬ 
tion of \Vus\f+'^{0) = 0.2165(4) from Moulson [26]® and the preliminary 2015 ( 2-|-1-1-1) -flavor 
FLAG average for /+(0)^^ = 0.9704(24)(22) [12, 13]^° gives |K,| = 0.22310(74)thy(41)exp 
from K £2 decay. The determinations of |14s| from leptonic and semileptonic kaon decays 
are both quite precise (with the error from leptonic decay being about 20% smaller), but 
the central values differ by 2.2a. Finally, the combination of the ratio |14s|/|Uud| from 
leptonic decays [Eq. (26)] with \Vud\ from f3 decay implies an alternative determination of 
|Uus| = 0.2254(6) which agrees with the value from leptonic kaon decay, but disagrees with 
the iC^s-decay result at the 2.2a level. Collectively, these results indicate that that there 
is some tension between theoretical calculations and/or measurements of leptonic pion and 
kaon decays, semileptonic kaon decays, and superallowed /3-decay. Although this may be 
due to the presence of new physics, it is also important to revisit the quoted uncertainties 
on both the theoretical and experimental inputs. 

Finally, we combine the experimental measurements of fT,+ \Vud\i fK+\^us\ from leptonic 
pseudoscalar-meson decays in Eqs. (10) and (11) with determinations of the CKM elements 
from other decays or unitarity to infer “experimental” values for the decay constants. As¬ 
suming that there are no signihcant new-physics contributions to any of the input processes, 
the comparison of these results with theoretical calculations of the decay constants enables 
a test of lattice-QCD methods. Taking |14d| from superallowed /3-decay [100] leads to 

^jexp” _ i 30.50(1)(3)(13) MeV, (28) 

where the uncertainties are from the errors on T, \Vud\, and higher-order corrections, respec¬ 
tively. This agrees with the theoretical value / 7 r+ = 130.2(1.7) MeV in Eq. (16) obtained 
from an average of recent (2-1-1)-flavor lattice-QCD results [35, 37, 39]. We take the value 
iKisI = 0.22534(65) from the most recent global unitarity-triangle ht of the UTfit Collabo¬ 
ration [101] because there is tension between the values of \Vus\ obtained from leptonic and 
semileptonic kaon decays. This implies 

= 155.72(17)(45)(16) MeV (29) 

® This is an update of the 2010 Flavianet review [28] that includes new measurements of the Ks lifetime [95, 
96], Re(e'/e) [96], and B{K^ —>• 7r=*=7r+7r“) [27]. The latter measurement is the primary source of the 
reduced error on via the constraint that the sum of all branching ratios must equal unity. 

This result comes from the calculation of FNAL/MILC in Ref. [97]. For comparison, the 2015 preliminary 
(2+l)-flavor FLAG average based on the calculations of FNAL/MILC [98] and RBC/UKQCD [99] is 
f+{0)^^ = 0.9677(37) . 
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where the uncertainties are from the errors on F, \Vus\^ and higher-order corrections, respec¬ 
tively. This agrees with the theoretical value /x+ = 155.6(0.4) MeV in Eq. (16) obtained 
from an average of recent three and four-flavor lattice-QCD results [31-33, 35, 37, 39]. 


B. \Vcd\, 114 sI, and status of second-row unitarity 

Using the average values for \Vcd\fD+ and \Vcs\f£)+ from Tables II and III, and for /£>+ 
and f^+ from Eq. (19), we obtain the following determinations of the CKM matrix elements 
\Vcd\ and | 14 s|, and from leptonic decays within the standard model: 

\VJ = 0.217(5)(1) and |14s| = 1.007(16)(5), (30) 

where the errors are from experiment and theory, respectively, and are currently limited by 
the measured uncertainties on the decay rates. The central value of |14s| is greater than 
one, but is compatible with unity within the error. The above results for \Vcd\ and |14s| 
do not include higher-order electroweak and hadronic corrections to the rate, in analogy to 
Eq. (3). These corrections have not been computed for Zl[^^-meson leptonic decays, but are 
estimated to be about to be about 1-2% for charged pion and kaon decays (see Sec. II A). 
Now that the uncertainties on 114^1 and |14s| from leptonic decays are at this level, we hope 
that the needed theoretical calculations will be undertaken. 

The CKM elements \Vcd\ and |I4s| can also be obtained from semileptonic —>■ 

and —)■ decays, respectively. Here experimental measurements determine the 

product of the form factor times the CKM element, and theory provides the value for the 
form factor at zero four-momentum transfer between the initial meson and the hnal 
pion or kaon. We combine the latest experimental averages for f^'"{0)\Vcd\ = 0.1425(19) 
and /+'"^(0)|I4s| = 0.728(5) from the Heavy Flavor Averaging Group (HFAG) [102] with 
the zero-momentum-transfer form factors /+’^(0) = 0.666(29) and /+'’^(0) = 0.747(19) 
calculated in (2-1-1)-flavor lattice QCD by the HPQCD Collaboration [103, 104] to obtain 
\Vcd\ = 0.2140(97) and |I4s| = 0.9746(257) from semileptonic Zl( 5 )-meson decays. The values 
of \Vcd\ from leptonic and semileptonic decays agree, while those for |I4s| are compatible 
at the l.lcT level. The determinations of |I4d| and |I4s| from leptonic decays in Eq. (30), 
however, are 2.0x and 1.6x more precise than those from semileptonic decays, respectively. 

The results for |I4d| and |I4s| from Eq. (30) enable a test of the unitarity of the second 
row of the CKM matrix. We obtain 

iKrfP + IK.P + iKfeP - 1 = 0.064(36), (31) 

which is in slight tension with three-generation unitarity at the 2 ct level. Because the 
contribution to Eq. (31) from \Vch\ is so small, we obtain the same result taking |I4bP“*' x 
10^ = 42.21(78) from inclusive B —)■ Xciv decay [105] or x 10^ = 39.04(75) from 

exclusive B —)■ D*£v decay at zero recoil [106]. 

We can also combine the experimental measurements of fD+\^cd\ = 45.91(1.05) MeV and 
Jd+I^csI = 250.9(4.0) MeV from leptonic pseudoscalar-meson decays from Tables H and HI 
with determinations of |I4d| and |I4s| from CKM unitarity to infer “experimental” values 
for the decay constants within the standard model. For this purpose, we obtain the values 
of \ Vcd\ and |I4s| by relating them to other CKM elements using the Wolfenstein parameter- 
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ization [107]. We take \Vcd\ to equal the value of \Vus\ minus the leading correction [108]: 


iWdI = IK, 
= |K 




-1 + (1 - 2p) 


IK 


cb \ 


+ p' 




\ 1/2 


(32) 

(33) 


Using |Ks| = 0.2255(3)(6)(3) from leptonic kaon decay [Eq. (26)], inclusive iKftl as above, 
and {p^p) = (0.136(24), 0.361(14)) from CKM unitarity [101] |Kd| =0.2254(7). We take 
IK^I = |Kd| - |KbP/2 [108], using |Kd| = 0.97417(21) from {3 decay [94], giving |Ks| = 
0.9733(2). Given these choices, we hnd 

/“K^” = 203.7(4.7)(0.6) MeV and /“+p” = 257.8(4.1)(0.1) MeV, (34) 


where the uncertainties are from the errors on E and |Ks| (or iKdl), respectively. These 
disagree with the theoretical values /£>+ = 211.9(1.1) MeV and f^+ = 249.0(1.2) MeV 
in Eq. (19) obtained from averaging recently published three and four-flavor lattice-QCD 
results at the 1.7 ct and 2.0a levels, respectively. The signihcances of the tensions are sensitive, 
however, to the choices made for iKsI and iKdl- Thus resolving the inconsistencies between 
determinations of elements of the hrst row of the CKM matrix discussed previously in 
Sec. V A may also reduce the mild tensions observed here. 


C. 114bI and other applications 

Using the average value for \Vub\fB+ from Table V, and for fB+ from Eq. (24), we obtain 
the following determination of the CKM matrix element |Kfe| from leptonic decays within 
the standard model: 

iKbl =4.12(37)(9) X 10-^ (35) 

where the errors are from experiment and theory, respectively. We note, however, that 
decays involving the third generation of quarks and leptons may be particularly sensitive to 
new physics associated with electroweak symmetry breaking due to their larger masses [4, 6] , 
so Eq. (35) is more likely to be influenced by new physics than the determinations of the 
elements of the hrst and second rows of the CKM matrix in the previous sections. 

The CKM element |Kb| can also be obtained from semileptonic S-meson decays. Over the 
past several years there has remained a persistent 2-3 ct tension between the determinations of 
I Kb I from exclusive B —)■ Trf'i/ decay and from inclusive B —)■ Xuf'z/ decay, where denotes 
all hadrons which contain a constituent up quark [3, 102, 109-111]. The currently most 
precise determination of iKbl*^^*^^ = 3.72(16) x 10“^ is obtained from a joint z-ht of the vector 
and scalar form factors and /(f’^(g^) calculated in (2-|-l)-havor lattice QCD by the 

FNAL/MILC Collaboration [112] and experimental measurements of the differential decay 
rate from BaBar [113, 114] and Belle [115, 116]. On the other hand, the most recent PDG 
average of inclusive determinations obtained using the theoretical frameworks in Refs. [117- 
119] is |Kb|™‘^’ = 4.49(16) (I}®) ^ 10”^ [^20]. The result for |Kb| from leptonic B ^ tv decay 
in Eq. (35) is compatible with determinations from both exclusive and inclusive semileptonic 
R-meson decays. 

The CKM element |Kb| can now also be obtained from semileptonic Ab decays. Specif¬ 
ically, the recent LHCb measurement of the ratio of decay rates for Ab —)■ p^v over Ab —)■ 
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Ac£z/ [121], when combined with the ratio of form factors from (2+l)-flavor lattice QCD [122], 
enables the hrst determination of the ratio of CKM elements |Vufe|/|V[;fe| = 0.083(4)(4) from 
baryonic decay. Taking = 42.21(78) x 10“^ [105] for the denominator,^^ we obtain 

iKiftl = 3.50(17)(17) (6) X 10“^ from exclusive A;, semileptonic decays, where the errors are 
from experiment, the form factors, and iKfel, respectively. The result for \Vub\ from leptonic 
B ^ TV decay in Eq. (35) is 1.4cr higher than the determination from 6-baryon decays. 

Given these results, the puzzle still stands, and the determination from leptonic B^- 
meson decay is not yet sufficiently precise to weigh in on the discrepancy. New and improved 
experimental measurements and theoretical calculations of other 6 —)■ u flavor-changing 
processes, however, are providing additional information and sharpening the picture of the 
various tensions. Further, the error on \Vub\ from B ^ tv decay will shrink once improved 
rate measurements from the Belle II experiment are available. 

Finally, we can combine the experimental measurement of |I4ib|/n+ from leptonic B^- 
meson decays in Table V with a determination of the CKM element \Vub\ from elsewhere to 
infer an “experimental” values for /s+ within the standard model. This, of course, assumes 
that there are no signihcant new-physics contributions to B^ — tv, which may turn out 
not to be the case. Further, one does not know a priori what value to take for \Vub\ given 
the inconsistencies between the various determinations discussed above. We therefore take 
the PDG weighted average of the determinations from inclusive and exclusive semileptonic 
B-meson decays |= 4.09(39) x 10“^ [120], where the error has been rescaled by 
the ^Jx^= 2.6 to account for the disagreement. Using this result we obtain 

4+P” = 188(17)(18) MeV, (36) 

where the uncertainties are from the errors on T and \Vub\i respectively. This agrees within 
large uncertainties with the theoretical value /b+ = 187.1(4.2) MeV in Eq. (24) obtained 
from an average of recent three and four-flavor lattice-QCD results [69, 86, 88, 89]. 
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